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ABSTRACT 


The  flow  field  of  a  sonic  jet  exhausting  counter  to  a  supersonic 
airstream  has  been  investigated  theoretically  and  experimentally. 
The  total  shock-layer  thickness  was  calculated  using  a  blunt-body- 
type  analysis  and  the  known  properties  of  a  free-jet  expansion  from 
a  sonic  orifice.  The  total  shock-layer  thickness  and  position  of  the 
outer  shock  relative  to  the  orifice  are  shown  to  be  a  function  of  free- 
stream  Mach  number,  jet  reservoir  pressure,  free-stream  pitot 
pressure,  and  orifice  size.  The  predicted  inner  shock  position  is 
compared  with  previously  published  experimental  data,  and  the  pre¬ 
dicted  outer  shock  position  is  compared  with  data  obtained  in  a  low 
density  wind  tunnel.  Results  in  the  transitional  flow  regime  indicate 
that  the  outer  shock  to  orifice  position  distance  is  greater  than  the 
predicted  value. 
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SECTION  I 
INTRODUCTION 


The  flow  pattern  of  a  jet  exhausting  counter  to  a  supersonic  free 
stream  has  been  investigated  by  several  authors  from  the  standpoint  of 
application  to  thrust  vector  control  and  thermal  protection  of  the  re-entry 
vehicle  (Refs.  1,  2,  and  3).  More  recently  the  case  of  a  sonic  jet  ex¬ 
hausting  from  the  nose  of  a  sounding  rocket  at  high  altitudes  has  become 
important.  As  part  of  the  AFCRL  upper  atmosphere  research  program, 
nitric  oxide  (NO)  was  released  from  the  nose  of  sounding  rockets  in  the 
upper  atmosphere  (90-  to  140-km  altitude)  at  night  to  determine  the  atomic 
oxygen  (O2)  concentration.  A  series  of  tests  simulating  the  rocket  flights 

was  performed  at  AEDC  to  determine  the  effective  rate  constant  (K)  for  the 

K 

reaction,  NO  +  0  :^r  NO2  +  hr.  Results  of  these  tests  are  reported  in 
Refs.  4,  5,  and  6.  This  report  presents  an  inviscid  analysis  of  the  flow 
field  of  a  sonic  jet  exhausting  counter  to  a  supersonic  free  stream  and 
discusses  the  effects  of  flow  rarefaction  on  the  flow  field.  Previous  ex¬ 
perimental  results  (Refs.  1,  2,  and  3)  are  compared  with  the  present 
theory,  and  data  obtained  from  a  low  density  wind  tunnel  are  presented. 

The  analysis  neglects  chemical  reactions  in  the  shock  layer. 


SECTION  II 

FREE-JET  FLOW  FIELD 


The  flow  field  of  the  sonic  jet  expanding  into  a  vacuum  has  been  in¬ 
vestigated  by  Ashkenas  and  Sherman  (Ref.  7),  who  obtained  fitting  formulas 
from  a  method  of  characteristics  solution  of  the  flow.  From  experimental 
pitot  probe  data,  Ashkenas  and  Sherman  obtained  the  following  formula 
along  the  centerline  for  2  <X/D  <  90: 

p'  /xy2-07  7 

— L=0.64°^-j  fory  =  T  (1) 

X  =  distance  along  centerline  from  orifice  exit, 

D  =  diameter  of  orifice, 

P c,'j  =  pitot  pressure  behind  a  normal  shock,  and 

Poj  =  jet  reservoir  pressure. 

Ashkenas  and  Sherman  also  compared  the  results  obtained  by  using  a 
sonic  nozzle  and  a  thin  plate  orifice  and  found  that  for  the  round  orifice. 
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the  effects  of  entrance  geometry  on  the  flow  field  were  negligible  for 
Ref)*  >  600,  where  Ref)*  equals  Reynolds  number  based  on  the  orifice 
diameter.  However,  viscous  effects  began  reducing  the  effective 
orifice  diameter  in  the  nozzle  orifice  at  Re^*  <  600.  This  point  will 
be  discussed  again  later  with  the  present  experimental  results  in 
Section  6.  2. 


SECTION  Hi 

FLOW  FIELD  CONFIGURATION 


A  detailed  sketch  of  the  inviscid  flow  field  is  shown  in  Fig.  1  and 
serves  to  define  symbols. 


3.1  general  structure 

The  general  configuration  has  been  investigated  (Refs.  1  and  2)  and 
consists  of  an  outer  shock  through  which  the  free  stream  is  decelerated, 
an  inner  shock  through  which  the  jet  flow  is  decelerated,  and  a  contact 
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surface  where  the  two  flows  meet  and  the  pressure  of  the  two  flows  at 
the  contact  surface  must  equilibrate.  The  shock  pattern  is  in  the  form 
of  an  approximate  spherical  shell  with  the  flow  between  the  two  shocks 
expanding  to  the  ambient  conditions  at  a  'shoulder  point"  located  where 
the  jet  boundary  meets  the  inner  shock. 


3.2  CONTACT  SURFACE  SHAPE 


It  can  be  shown  that  the  contact  surface  is  approximately  spherical 
by  comparing  the  pressure  distribution  on  the  free-stream  side  of  the 
contact  surface  with  the  pressure  on  the  jet  side.  The  pressure  distribu¬ 
tion  on  the  free-stream  side  is  obtained  from  a  modified  Newtonian 
analysis  with  the  contact  surface  assumed  a  solid  sphere,  or 


-^7-  =  cos2  $  +  -~r-  sin!  $ 

*900  *  9  i-Vi 


(2) 


where 


P  =  pressure  at  contact  surface, 

p2«=  =  pressure  at  stagnation  point  on  contact  surface, 
and 

-  free-stream  static  pressure 

which  is  plotted  in  Fig.  2.  The  pressure  distribution  on  the  jet  side  is 
ohtained  by  using  the  formula  derived  by  Sherman  {Ref.  8)  for  the  density 
ratio  as  a  function  of  $ , 

=  ■  ■■■-  —  cos2  (0.945$)  for  y  -  —  (3) 

poj  R1  5 

Pil  -  static  density  in  the  free  jet, 

Poj  *  jet  reservoir  density,  and 
H  =  distance  from  orifice  exit 

from  which  the  free-jet  Mach  number  and  pressure  ratio  across  the  inner 
shock  (assumed  normal)  can  be  calculated  from  isentropic  relations.  The 
calculated  distribution  is  compared  with  the  modified  Newtonian  result  in 
Fig.  2.  The  pressure  distributions  agree  very  closely,  confirming  the 
approximate  spherical  nature  of  the  contact  surface. 
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3.3  INNER  SHOCK  LOCATION 

Previous  experimental  data  for  continuum  air  for  the  inner  shock 
location  on  the  jet  axis  are  summarized  in  Fig.  3  along  with  Charwat’s 
(Ref.  3)  theory  for  the  inner  shock  location.  The  experimental  data 
(Refs.  2  and  3)  were  obtained  from  schlieren  photographs  at  Mro  =2.9 
to  7.  1. 


Fig.  3  Inner  Shock  Location 


f 
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The  inner  shock  location  can  also  be  found  from  Eq.  (1)  by  assuming 
that  the  free-jet  pitot  pressure  behind  a  normal  shock  equates  with  the 
free-stream  pitot  pressure  along  the  jet  centerline;  i.  e.  , 

P '  ■  -  P" 

*  OJ  1  QOO 


X 

D 


D 


0.806 


(4) 


which  is  plotted  in  Fig.  3.  Equation  (4)  does  give  somewhat  better 
agreement  with  the  experimental  shock  locations  and  will  be  used  in 
further  derivations. 


SECTION  IV 

SHOCK-LAYER  THICKNESS 

4.1  INVISCID  SOLUTION  FOR  THE  FLOW  FIELD  BETWEEN  THE  CONTACT  SURFACE 
AND  THE  INNER  SHOCK 

The  flow  field  between  the  inner  shock  and  the  contact  surface  is 
analogous  to  a  blunt-body  flow  field  where  the  shock  is  detached  from 
the  surface.  Upon  analysis  of  the  known  parameters,  it  is  found  that 
a  hypersonic  blunt-body  analysis,  similar  to  Li  and  Geiger's  constant 
density  solution  (Ref.  9)  can  be  adapted  to  the  present  flow  field  to  solve 
for  the  shock-layer  thickness  at  the  stagnation  point.  A  correlation 
between  free-stream  conditions,  jet-stream  conditions,  jet  total  condi¬ 
tions,  and  shock  detachment  distance,  Aj,  can  be  made.  Previous 
experimental  data  on  blunt  bodies  in  hypersonic  flow  (Ref.  9)  have  shown 
Li  and  Geiger's  solution  to  be  accurate  for  hypersonic  Mach  numbers, 
which  is  the  case  for  the  jet  flow  fields  considered  here. 

The  following  derivation  and  discussion  review,  in  part,  Li  and 
Geiger's  derivation  and  apply  boundary  conditions  and  assumptions 
appropriate  for  the  present  flow  field.  Symbols  are  defined  in  Fig.  1. 

The  basic  equations  of  motion  governing  the  axisymmetric  flow  of  an 
inviscid  fluid  are: 


1 

v  da  + 

K 

(1  +  Ky) 

u  + 

dy 

(1  +  Ky) 

+  1  -  i  =  0 

(1  -t-  Ky)  p  <?x 


(5) 
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1 

dv  dv 

K 

(1  +  Ky) 

d*  "  V  dy 

( 1  +  Ky) 

+ 


1  dP  n 

—  3—  «  u 

p  Oy 


1 

O  +  Ky  > 


J-  ipv)  + 
dy 


K 

(1  +  Ky) 


pv 


+ 


J_  I  P” 
r  |[jl  +  Ky) 


dr 

d7 


+ 


0 


(6) 


(7) 


where  x  and  y  are  distances  measured  along,  and  normal  to,  the  contact 
surface,  K  is  the  curvature  of  the  contact  surface,  or  K  =  1/RC. 


Certain  simplifications  can  be  applied  to  the  governing  equations; 

1.  Since  Ai  is  small  in  comparison  to  Rc  and  Ky  =  0(A4/Rc),  then 
Ky  «  1. 

2.  Since  the  velocities  in  the  shock  layer  at  the  stagnation  point 
(u,  v)  are  small,  then  terms  Kuv  and  Ku2  can  be  neglected. 

3.  In  the  vicinity  of  the  stagnation  point,  r  «  x.  Therefore,  for 
small  r,  r  can  be  replaced  by  x. 

4.  As  M}j—  »,  M2j— *-(7-  1  >/ 2-y.  Therefore,  the  flow  field  near 

the  stagnation  point  can  be  regarded  as  essentially  incompressible, 
or  P  -  P 2j  =  P lj/e  • 


Applying  these  simplifications  to  Eqs.  (5),  (6),  and  (7),  the  governing 
equations  become; 


u 


da 

dx 


a 


d\ 

dx 


d(ux)  + 
dx 


du  €  dP 

dy  pVi  dx 

(8) 

v  —  -  —  f  dP 
dy  pls  dx.^ 

(9) 

d(vx)  Q 

dy 

{10) 

The  appropriate  boundary  conditions  are: 

1.  Stagnation  conditions  on  the  contact  surface,  y  =  0,  x  -  0,  then 
u  =  v  =  0. 
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2.  Condition  of  tangency  along  the  contact  surface  is  y  =  0,  v  =  0. 

3.  Conditions  along  the  inner  shock  are  y  =  Aj,  u  =  0, 

v2j  ‘ 

4.  Conditions  along  the  contact  surface  were  defined  in  the  previous 
section  by  matching  the  flow  in  the  inner  shock  layer  with  the 
flow  in  the  outer  shock  layer.  It  was  shown  that  the  pressure 
distribution  was  equivalent  tc  the  modified  Newtonian  distribution. 


or 


cp  =  cp0  cos*  ^ 

p  p, 

- —  =  cos2  $  +  -  sin2  <t> 

A  i  •  n  i 


(ID 


Tne  incompressible  Bernoulli  equation  in  the  stagnation  region  is: 


°]  ~  P  +  9  Poj  1)2 


or 


\ 


2p; 


P  o3 


1  ~ 


P' 

roj 


(12) 


(13) 


Substituting  Eq.  (11)  into  Eq.  (13), 


u  =  s  in  4’ . 


2V 


~pi V 


1  - 


Ij 


®J 


(14) 


If  4>  is  small,  then  sin  «  <f>  and  x  =  R j  <J> .  Therefore, 


/  2P  ' 
/  ^  o  1 


“  -  hT 


1  - 


p' 

°) 


(15) 


4.2  SOLUTION  FOR  THE  INNER  SHOCK  DETACHMENT  DISTANCE 

The  P  is  first  eliminated  from  Eqs.  (8)  and  (9)  by  differentiating 
Eq.  (8)  with  9/9y  and  Eq.  (9)  with  9/dx,  or 

dn  dj  +  u  +  dv  da  +  v  £» 

dy  dx  dxdy  dy  dy  dy2 


f  d2P 

p  ij  dxdy  (16) 
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du  dv  +  u  d2v  +  dy.  _3v  +  v  d2v 
dx  dx  dxJ  dx  dy  dxdy 


_  <yp 

P U  dxdy 

By  subtracting,  one  obtains: 


dj  du  t  u  d2u  +  .dy  jdu  +  v  d2u  ___  du.  _dv 

dy  dx  dxdy  dy  dy  dy2  dx  dx 


d  v  _  _dv  dy  _  dav 
(3xa  <9x  dy  dxdy 


Combining  terms  in  Eq.  (12): 


<9*u  d2v 

d2*  _  d*v 

dxdy  dx 2 

+  V 

dy 2  dxdy 

du  _ 

dv 

du  _  dy 

dy 

dx 

+  ~ 
dx 

dy  dx 

(17) 


(18) 


(19) 


The  form  of  u  and  v  must  be  chosen,  and  the  following  argument  is 
presented  by  Li  and  Geiger  (Ref.  9).  For  small  x,  u(x,  y)  may  be 
expressed  as  an  ascending  series  of  the  form, 

u(x,y)  -  aQ  +  a [X  +  a2x2  +  ,  .  .  (20) 

where  the  a's  are  functions  of  y.  The  symmetry  of  the  problem  permits 
only  terms  in  odd  powers  of  x,  and  for  small  x's  third  and  higher  order 
terms  are  unimportant.  Therefore,  it  can  be  assumed  that 


U  =  X  fly)  (21) 

where  f'  -  df/dy,  or  a  linear  variation  of  u  with  respect  to  x  near  the 
stagnation  point.  Substituting  Eq,  (21)  into  the  continuity  Eq.  (10),  the 
form  of  v  is: 


v  =  -2  f(y) 


(22) 


Equations  (19),  (21),  and  (22)  are  now  combined  in  the  following 
manner: 

xf'  0)  -  2f  (xf"'  -  0)  -  2f '  (xf"  -  0)  +  r  (xf"  -  0)  =  0 

where  f°  =  d^f/dy*and  f"'  =  d^f/dy^,  or 
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-2xll'"  =  0 

and  since  f  4  0,  then 


l'"  =  0 


(24) 

(25) 


Equation  (25)  is  a  simple  differential  equation  which  can  be  solved 
to  yield; 

f(y)  =  C0  +  C,y  +  C2y2 - j  (26) 

where  the  constants  of  integration  must  be  determined  from  the  boundary 
conditions. 

The  condition  of  tangency  (y  =  0,  v  =  0)  yields  C0  =  0. 

The  constants  Cj  and  C2  are  determined  from  the  conditions  at  the 
inner  shock. 

As  x— 0 


y  =  Aj 


or 

f(Aj)  -  CtAj  +  C2AJj 

but 

f(A.)  .  + 

and,  since 

Pi\ 

v  c  ~V2J  =  ~  PIT  v*i  =  "£v‘i 

or 

C  jAj  +  CjAJj  =  +  ^  v,j 

From  Eq.  (21), 

fly)  =  C.  +  2Csy  =  41- 


and,  as  x— -0,  y— - Aj;  therefore, 

C,  +  2CJAJ  =  x^0  (-£-) 


but  from  Eq.  (15), 


(27) 


(28) 


(29) 
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Substituting  Eq.  (28)  into  Eq,  (29): 


C,  +  2CjA| 


(30) 


Equations  (27)  and  (30)  can  now  be  solved  simultaneously  for  Cj  and  Cg 
in  terms  of  Aj,  Rj,  P'Qj,  p 'Qy  P,  e  ,  and  v^.  Or, 


C,  -  + 


(31) 


Ca 


«V, 


2'Vj 


1  I  2 P'oj 


(32) 


To  solve  for  Aj,  9P/8x  is  found  from  Eq.  (11)  and  substituted  into 
Eq.  (8)  along  with  Eqs.  (21)  and  (22).  Then  can  be  found  in  terms 
of  Rj,  P'Qj,  p’0-,  and  P^.  From  Eq.  (11), 


and 


or 


(33) 


(34) 


as  x — 0,  (sin  4>)/x — -1/Rj,  cos  <t> — 1 ►!,  and  (9<I>/8x) — *-l/Rj.  Therefore, 


(35) 


Substituting  Eqs.  (26)  and  (32)  into  (35)  and  solving  for  Cj: 


(36) 
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Substituting  Eq.  (36)  into  Eq.  (31)  and  choosing  the  (+)  sign  in  Eq.  (36), 


but 


and  since 


Therefore, 


^ j m  ax 


1  for  M,j  1 


and  also. 


and  since 
1 


Then, 


1  2  y 

vlj 

\  y  - 1 

l  m  / 

* 

mP'of 

mP  0 , 

P  Oj 

“  A  , 

RT'0J 

RT0) 

M 

|  -  1,  and  ( 

=  7* 

P'oj 

=  0.156Rj 

for  y  =■  ] 

since  TOJ  r  T'0j 


(37) 


(38) 


Equation  (4)  relates  to  the  ratio  of  jet  reservoir  pressure  to  free- 

stream  pitot  pressure.  Therefore,  substituting  Eq.  (4)  into  Eq.  (38) 
yields: 


Aj  =  0.1257D 


or 


Rc  -  Rj  -  Aj  =  0.932D 


(39) 


(40) 
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Equation  (40)  shows  that  the  contact  surface  location  is  related  to 
the  orifice  size,  D,  and  to  the  ratio  of  the  jet  reservoir  pressure/ 
free-stream  pitot  pressure  (P^/P^).  Note  that  Eq.  (40)  produces 
accurate  predictions  for  hypersonic  free-jet  Mach  numbers,  i.  e.  , 

Mjj  >  5,  From  Ashkenas  and  Sherman  (Ref.  7)  Mjj  =  5  occurs  at 

approximately  x/D  =  3.1  or  from  Eq.  (1),  P0j/P  ^  =  16.  3  which  should 
be  considered  the  lower  pressure  ratio  limit ’’of  applicability  of  Eq.  (40). 


4.3  DETACHMENT  DISTANCE  OF  OUTER  SHOCK  FROM  CONTACT  SURFACE 

The  distance  from  the  contact  surface  position  to  the  outer  shock  is 
determined  by  assuming  the  contact  surface  to  be  an  equivalent  spherical 
body  and  calculating  an  equivalent  shock  detachment  distance  using  avail¬ 
able  numerically  calculated  results,  Eq.  (10).  For  hypersonic  flow  the 
constant  density  assumption  for  the  shock  layer  is  an  accurate  approxima¬ 
tion  since  the  Mach  numbers  in  the  shock  layer  are  low  (from  P\! f>2 - *- 

small),  and  the  density  is  essentially  constant  from  the  shock  to  the  body. 
However,  for  supersonic  flow  the  density  ratio  is  no  longer  small,  and 
velocities  in  the  shock  layer  are  relatively  high,  producing  large  density 
gradients  between  the  shock  and  body  (contact  surface).  Hence,  the 
numerical  calculations  of  Van  Dyke  (Ref.  10),  which  show  good  agreement 
with  experiment  in  the  low  Mach  number  range,  are  used  to  determine 
A«,  (Fig.  4). 


Fig.  4  Outer  Shock  Detachment  Distance  from  Contact  Surface 
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Therefore,  the  distance  from  the  orifice  exit  to  the  outer  shock  is: 


Re 


(41) 


=  0.932D 


(42) 


where  A^/Rc  is  found  for  M^,  from  Fig.  4. 


SECTION  V 
EXPERIMENTS 


5.1  TEST  CHAMBER 

The  tests  were  conducted  in  the  Aerospace  Research  Chamber  (8V), 
(ARC  8V)  at  AEDC.  The  pumping  system  for  the  ARC  8V  consists  of  a 
6-in.  oil  diffusion  pump  backed  by  a  140-cfm  mechanical  pump  which  is 
used  for  rough. pumping  of  the  chamber.  During  continuous  runs,  the  gas 
flow  is  cryogenically  pumped  on  liquid  nitrogen  (LN2)  (77°K)  and  gaseous 
helium  (GHe)  (20°K)  cryosurfaces  totaling  1140  ft2.  Run  times  from  10  to 
20  min  could  be  made  at  normal  operating  mass  flow  rates.  A  nominal 
Mach  number  3  open  test  section  nozzle  was  used  to  produce  the  supersonic 
free  stream  against  which  the  sonic  jet  exhausted  (Fig.  5). 


L^-Cooled  Nozzle 
M-  3. 2-3. 5 
d'  *  27  cm 

dexlt "  76  cm  1 


Tubular  Gas 

Heater 

(29Q-100CPK) 


Model  with  NO  Radja||y  positioned 
Supply  Une-7  Cryopane|s  GHe 
/  Cooled  (15-20PK) 


Plenum  Chamber 
(Heated) 

p  -  50-1000  pHg 
T0  *  290-100CPK 


S~L 


35-mm  Camera 


//Pressure! 

//Transduce 


. . .  .  ...  / 

•  o  •  jL 

©  a  e  a  ©jo  «  •  \\r\ 


LN2-Cooled  End 
Panel  and  Liner 
Cryopanels 


'  LN2-Cooled  Circular 
Duct  and  Cone 


Fig.  5  Top  View  Schematic  of  Low  Density  Wind  Tunnel  (ARC  8V) 
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A  mixture  of  approximately  80  percent  Ng  -  20  percent  O2  was  used  for 
the  supersonic  nozzle  flow.  Three  models  were  tested.  The  models 
were  designed  so  that  an  adequate  (5-  to  10-cm)  size  outer  shock  radius 
would  be  produced  over  the  range  of  pressure  ratios,  P„-j/P'  =  10  to 

*  **  Odd 

10’.  Test  conditions  are  listed  in  Table  I.  Nitric  oxide  was  used  in 
the  models  since  this  gas  was  being  used  in  a  simultaneous  investigation 
of  the  chemiluminescent  reaction  of  NO  +  O2  reported  in  Ref.  9. 

TABLE  I 

WIND  TUNNEL  CONDITIONS 


po®'  M  Hg 

P'cwMHg 

R  1 

98.  6 

25.  5 

290 

3.  27 

0.  54 

8.  9 

27 

198.  5 

47.  0 

290 

3.  38 

0.  29 

17.  2 

40 

298.  5 

67.  0 

290 

3.  44 

0.  202 

25.  0 

45 

379.  0 

81.  5 

290 

3.  49 

0.  165 

50 

104.  0 

29.  0 

700 

3.  19 

1.  60 

20 

301.  0 

70.  0 

700 

3.  40 

0.  64 

8.  0 

30 

5.2  NOZZLE  CALIBRATION 

The  nozzle  was  calibrated  with  a  pitot  probe  over  the  range  of 
reservoir  pressures  (P0cc)  and  total  temperature  (T0m)  used  for  the 
tests.  The  actual  Mach  number  and  other  free- stream  parameters  are 
listed  in  Table  I.  The  Mach  number  variation  is  caused  by  the  change 
in  nozzle  boundary-layer  thickness  as  the  nozzle  flow  density  changes 
with  reservoir  conditions.  The  usable  test  core  diameter  varied  from 
20  to  50  cm  for  the  test  conditions  listed. 


5.3  MODELS 

Three  models  were  tested  (Fig.  6).  The  orifice  diameters  were 
0.  05,  0.  15,  and  0.  5  in. ,  respectively.  Reservoir  pressure  (PQj)  in 
the  models  was  measured  with  strain- gage  transducers  on  the  two 
smaller  models  (D  =  0.  05,  0.  15  in.  ,  PGj  =  5  to  220  mm  Hg)  and  with  a 
capacitance  transducer  on  the  largest  model  (D  =  0.5  in. ,  PQi  =  0.  5  to 
2.  2  mm  Hg). 
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Orifice 


Orifice 


c.  Model  3 

Fig.  6  Models  1,  2,  and  3 
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Nitric  oxide  was  supplied  to  the  models  from  standard  gas  bottles, 
necessary  valving,  and  flowmeters  outside  the  chamber.  To  prevent 
the  model  gas  supply  line  from  being  cooled  by  radiation  to  the  cold 
chamber  walls,  heaters  on  the  line  maintained  the  temperature  at  290°K. 

5.4  FLOW  VISUALIZATION 

Radio-frequency  (R-F)  excitation  (N2  afterglow)  was  used  to  visualize 
the  gross  structure  (i.  e.  ,  shocks)  of  the  flow  field.  Two  R-F  generators 
were  used.  The  R-F  signal  produced  by  the  two  generators  caused  the 
flow  field  to  be  illuminated,  and  the  outer  shock  could  be  clearly  distin¬ 
guished  for  most  of  the  flow  conditions.  One  antenna  was  attached  to  the 
insulated  copper  ring  at  the  nozzle  exit  plane.  The  other  antenna  was 
attached  to  the  model,  which  was  insulated  from  its  supports  by  a  nylon 
insert  (Fig.  7).  The  chamber  was  the  ground  for  both  generators. 


Pressure  Transducer 


,  10  llL. 

Scale 


Fig.  7  Model  Installation  ond  Flow  Visualization 

5.5  PHOTOGRAPHIC  TECHNIQUE  AND  ANALYSIS 

Photographs  of  the  model  and  illuminated  flow  field  were  made  through 
a  chamber  port  with  a  35-mm  camera  with  a  f/1.  8  lens  using  black  and 
white  Kodak®  Plus-X  film.  Figure  8  shows  a  typical  photograph  of  the  flow 
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configuration.  Microdensitometer  traces  of  the  film  negatives  (for 
example.  Fig.  9)  were  made  to  determine  the  relative  positions  of  the 
shock  and  the  orifice  exit  along  the  centerline  of  the  model.  The  model 
diameter  was  used  as  a  reference  dimension  on  the  negatives.  A 
definite  shock  position  was  difficult  for  the  T0m  =  700°K  runs  because  of 
the  gradual  film  negative  density  change  caused  by  the  thick  shocks  at 
these  low  densities. 


Fig.  9  Microdonsitometer  Trac« 
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SECTION  VI 

RESULTS  AND  DISCUSSION 

6.1  EXPERIMENTAL  RESULTS 

The  location  of  the  outer  shock  relative  to  the  orifice  position  is 
nondimensionalized  by  the  respective  orifice  diameter  and  is  plotted 
versus  P0]/Poa>  in  Pig'  10  for  all  the  runs.  Since  the  Mach  number 
range  (M^)  in  the  tests  was  relatively  low,  all  the  data  have  been 
plotted  together  so  that  the  overall  trend  in  the  data  can  be  seen.  Equa¬ 
tions  (4),  (40),  and  (42)  are  also  indicated  along  with  the  predicted  outer 
shock  location  using  Fig.  4  for  the  Mach  number  range  in  the  experiments. 
The  results  show  (Re/D)ex  >  (Re/D)^eo  1°  to  20  percent  fqr  all  condi¬ 
tions  tested.  ^ 

G 

Q  Oi/1 


Fig.  10  Experiment  I  Outer  Shack  Positions 

Neither  the  contact  surface  nor  the  inner  shock  could  be  seen  during 
the  runs  or  on  the  photographs.  Since  a  thorough  investigation  of  the 
excitation  and  emission  mechanisms  of  the  glow  is  beyond  the  scope  of 
this  report,  it  can  only  be  suggested  that  flow  inside  the  free-jet  boundary 
shock  was  not  excited  sufficiently  to  produce  a  distinct  glow  intensity 
change  at  the  inner  shock  position. 
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Since  neither  the  contact  surface  nor  the  inner  shock  could  be  seen, 
it  is  difficult  to  ascertain  a  priori  if  the  discrepancy  is  caused  by  error 
in  the  theoretical  prediction  or  by  rarefaction  effects  on  the  flow  field. 

Note  that  Rj/D  from  Eq.  (4)  is  approximately  5  percent  less  than  previous 
experimental  data  in  Fig.  3,  which  could  account  for  some  of  the  discrep¬ 
ancy  here.  By  examining  other  works  (Refs.  11,  12,  and  13)  on  low  density 
gas  dynamic  effects,  the  present  data  trend  can  be  seen  to  be  logical  for  the 
Reynolds  number  range  of  the  experiment. 


6.2  FLOW  RAREFACTION  EFFECTS 

For  the  case  where  the  molecular  mean  free  path  is  small  in  com¬ 
parison  to  a  characteristic  dimension,  the  gas  may  he  considered  a 
continuum,  and  gross  changes  in  the  flow  field  (i.  e.  ,  shock  waves)  occur 
over  a  very  short  distance.  However,  if  the  mean  free  path  is  of  the  same 
order-of-magnitude  as  the  characteristic  dimension,  gross  changes  in  the 
flow  field  may  occur  over  distances  comparable  to  the  characteristic 
dimension.  The  rarefied  flow  field  investigated  herein  differs  from  the 
idealized  inviscid  case  considered  in  Section  IV  in  that  the  inner  and  outer 
shocks  can  no  longer  be  considered  discontinuities,  and  the  shock  thick¬ 
nesses  become  significant  in  comparison  with  the  total  shock-layer 
thickness. 

Again  considering  the  similarity  of  the  present  flow  field  and  a  blunt- 
body  flow  field,  it  can  be  shown  that  the  effects  of  flow  rarefaction  in  the 
present  flow  field  would  be  similar  to  low  density  effects  in  blunt-body 
flow.  The  viscous  shock  layer  on  blunt  bodies  in  rarefied  flow  has  been 
discussed  in  Refs.  11  and  12.  As  the  Reynolds  number,  based  on  sphere 
radius  and  viscosity  downstream  of  the  normal  shock,  decreases  below 
=  1000,  the  stagnation  point  shock  detachment  distance  at  first  decreased 
slightly  and  then  increased  rapidly  below  Re  =  1000.  Bailey  and  Sims 
(Ref.  13)  confirmed  this  trend  for  spheres  and  flat-faced  bodies  in  low 
density,  hypersonic,  argon  (Ar)  flow.  These  trends  in  the  shock  detach¬ 
ment  distance  for  spheres  in  the  transition  regime  can  be  applied  to  the 
flow  field  considered  here  if  the  boundary  conditions  are  similar  for  both 
cases. 

Note  that  the  boundary  conditions  for  the  viscous  shock-layer  stagna¬ 
tion  point  solution  are  basically  the  same  as  the  boundary  conditions  for 
the  present  flow  field;  i.  e. ,  the  normal  velocity  component  at  the  contact 
surface  (wall)  =0,  and  the  conditions  behind  the  shock  are  given  by  the 
Rankine-Hugoniot  relations  or  a  Navier-Stokes  analysis.  The  Rankir.e- 
Hugoniot  relations  give  the  conditions  behind  the  normal  shock  for  the 
thin  shock  case.  For  the  thick  shock  case,  the  flow  field  from  the  free 
stream  (or  jet  flow)  to  the  contact  surface  must  be  completely  coupled  with 
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a  system  of  equations  such  as  the  Navier  Stokes  equations  (Ref.  12). 
Figure  11  shows  the  low  density  flow  field. 


Table  II  shows  the  Reynolds  numbers  for  the  jet  and  free-stream 
flow  fields  based  on  the  calculated  inviscid  shock  detachment  distance 
(Eq.  (42)  and  Fig.  4)  and  viscosity  behind  the  normal  shock, 


Re 


and 


He  >  _  p'll]Ai 

jA‘  ■ 


Remembering  that  for  hypersonic  flow  (Ref.  14), 

A _  , 


R 


then 


e  qu  Iv 


D-xi 

CAX  = - - - -  ~ 


^  e  qu  iv 


(43) 
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and 


RejAj 


^j^jRequiv  t 


-  Re 


R 


e 

c  q\i  tv 


(44) 


the  low  Reynolds  number  effects  on  the  blunt-body  shock  detachment 
distance  can  be  applied.  Therefore,  low  Reynolds  number  effects 
would  be  possible  in  the  present  case  for 


ReooAoa  =  RejAj  <  170,  y  =  1.4  (45) 

with  rapid  increases  in  A  for 

ReooA«  =  ReA  <  17,  y  “  1.4 

The  Reynolds  numbers  listed  in  Table  II  are  then  in  the  range  where 
increases  in  the  Shock  detachment  distances  would  be  expected. 
Rarefaction  of  the  jet  and/or  the  free-stream  flow  would  cause  the 
observed  increase  in  Re/D.  Normally,  one  would  expect  rarefaction  of 
both  flow  fields  simultaneously  unless  gases  of  widely  differing  properties 
are  used. 


TABLE  II 

MODEL  CONDITIONS  AND  FLOW  REYNOLDS  NUMBERS 


Model 

PQj,  mm  Hg 

O 

■'“3 

0 

ReD* 

ReJ^i 

A 

1 

44  to  223 

290 

1100  to  5580 

40.  8  to  74.  2 

1.113  to  10.  58 

2 

5  to  22 

290 

375  to  1650 

33.  1  to  40.  3 

1.  14  to  10.  40 

3 

0. 5  to  2.  2 

290 

125  to  550 

34.  2  to  39. 2 

1.  19  to  11. 58 

The  assumption  of  a  distinct  contact  surface,  as  was  done  for  the 
inviscid  solution,  would  not  be  valid  as  the  density  in  the  shock  layer  is 
reduced.  Considering  the  flow  in  the  shock  layer  on  a  microscopic  scale, 
the  migration  of  one  species  of  gas  across  the  "contact  surface"  would 
increase  as  the  density  is  reduced  and  mean  free  path  increases.  As  both 
shocks  become  very  thick,  the  shock  layer  would  become  fully  merged. 
The  analysis  would  then  necessarily  be  completely  coupled  to  determine 
species  concentration  profiles  in  the  shock  layer. 

Viscous  effects  at  the  orifice  should  also  be  mentioned.  The 
Reynolds  number  ranges,  based  on  sonic  conditions  and  the  orifice  diam¬ 
eter  (Reg)*  =  p*U*D*/p*),  are  shown  in  Table  II.  Results  by  Ashkenas 
and  Sherman  (Ref.  7)  and  Smetana,  et  al.  (Ref.  15)  for  a  similar  type  of 
nozzle  indicated  viscous  influences  on  the  "effective  '  orifice  diameter  to 
be  present  for  Re^*  <  600.  Based  on  this  estimate,  viscous  influences 
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on  the  effective  orifice  diameter  would  be  expected  for  Models  2  and  3. 
If  an  effective  orifice  diameter  has  been  used  in  Eqs.  (4),  (40),  and  (42) 
and  for  nondimensionalizing  the  experimental  data,  both  the  theoretical 
lines  and  the  data  points  would  be  shifted  equally,  and  the  trend  would 
remain. 


SECTION  VII 
CONCLUDING  REMARKS 


A  method  for  calculating  the  inviscid  flow  field  structure  of  a  sonic 
jet  exhausting  counter  to  a  supersonic  airstream  has  been  developed 
along  the  centerline.  The  predicted  inner  shock  to  orifice  exit  distance 
shows  good  agreement  with  previously  published  data  in  the  continuum 
flow  regime.  The  data  reported  herein  which  were  obtained  in  the  transi¬ 
tion  flow  regime  indicate  that  the  outer  shock  to  orifice  exit  distance  is 
greater  than  the  predicted  inviscid  value.  The  trend  appears  to  be  cor¬ 
rect  for  this  flow  regime. 

Several  recommendations  for  future  study  can  be  made  based  on  the 
present  experiment: 

1.  Extend  the  test  regime  to  lower  and  higher  densities  so  that 
continuum  regime  data  can  be  obtained  to  better  confirm  the 
theory  and  low  transition  regime  data  can  better  define  the 
trend  with  decreasing  Reynolds  numbers. 

2.  Develop  techniques  to  locate  the  inner  shock  position  for  all  flow 
regimes,  possibly  using  an  electron  beam. 

3.  Investigate  the  effects  of  mixing  at  the  contact  surface  using 
gases  with  widely  differing  thermodynamic  properties. 
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